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Abstract 

Evanescent random fields arise as a component of the 2-D Wold 
decomposition of homogenous random fields. Besides their theoretical 
importance, evanescent random fields have a number of practical ap- 
plications, such as in modeling the observed signal in the space time 
adaptive processing (STAP) of airborne radar data. In this paper we 
derive an expression for the rank of the low-rank covariance matrix of 
a finite dimension sample from an evanescent random field. It is shown 
that the rank of this covariance matrix is completely determined by 
the evanescent field spectral support parameters, alone. Thus, the 
problem of estimating the rank lends itself to a solution that avoids 
the need to estimate the rank from the sample covariance matrix. We 
show that this result can be immediately applied to considerably sim- 
plify the estimation of the rank of the interference covariance matrix 
in the STAP problem. 
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1 Introduction 



1.1 The Evanescent Random field 

The problem of linear prediction of stationary processes is a classic problem 
in time-series analysis. One of the most fundamental results in this field 
is the Wold decomposition [T3], that states that a regular one dimensional 
wide-sense stationary processes indexed by Z may be decomposed into two 
stationary and orthogonal components: the purely-indeterministic process 
(that produces the innovations) and the deterministic process. This decom- 
position can be equivalently reformulated using spectral notations: the spec- 
tral measure of the purely-indeterministic process is absolutely continuous 
with respect to the Lebesgue measure, and the spectral measure of the de- 
terministic process is singular. In other words, the spectral measures of the 
orthogonal components of Wold decomposition yield the Lebesgue decompo- 
sition of the spectral measure of the process. 

Homogenous random fields, (also called doubly stationary series), are 
the two-dimensional (indexed by Z 2 ) generalization of one-dimensional wide- 
sense stationary process. Unfortunately, unlike the one-dimensional case, in 
multiple dimensions there is no natural order definition and terms such as 
"past" and "future" are meaningless unless defined with respect to a specific 
order. Linear prediction of homogenous random fields was first rigorously 
formulated by Helson and Lowdenslager in [Hj. The problem of defining 
"past" and "future" on the two-dimensional lattice (i.e., Z 2 ) was resolved in 
[6] in terms of "half plane" total-ordering. The trivial example of a half-plane 
total order on Z 2 is a usual lexicographic order: (k, I) ^ (n, m) iff' k < n or 
(k = n and I < m). Lexicographic order can be considered as a linear order 
induced by Non-Symmetric (delimited by a broken straight line) Half Plane 
(NSHP), (see Figure |T]). 

Further analysis of the prediction problem led to a generalization of the 
Wold decomposition [7]. When we consider random processes indexed by a 
group we obtain a Wold decomposition with respect to any given total order 
on the group. When the group is not Z (like R or Z 2 ) the deterministic 
process can have as a direct summand a deterministic process of a special 
type, the evanescent process. In order to provide some intuition on the char- 
acteristics of the evanescent process we next state some basic definitions and 
present an example of an evanescent field defined with respect to a vertical 
total order, which is simply a lexicographic order on Z 2 : 
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Figure 1: Non-symmetric half plane 



A homogeneous random field {y(n,m)} is called regular with respect to 
the lexicographic order if for every (n,m), E[y(n,m) — y(n,m)] 2 = a 2 > 

where y(n, m) is the projection of y(n, m) on the c.l.m. {y(k, I) : k < n, I G 

Z}U{y(n,l) : I < m} , where c.l.m. denote a closed linear manifold. Thus, a 

regular homogeneous random field has a non-zero innovation at every lattice 
point. A homogeneous random field {z(n,m)} is called deterministic with 
respect to the lexicographic order if it can be perfectly linearly predicted 
from its past in mean-square sense, i.e., for every (n,m) we have z(n,m) G 

c.l.m. {z{k, I) : k < n, I G Z} U {z{n, I) : I < m} . 

Although a deterministic field can be perfectly predicted from its past 
with respect to lexicographic order, it may still posses a non-zero innovation 
when prediction is based on samples in previous columns only. We then 
say that the field {z(n, m)} has vertical column-to- column innovations if 
I(n,m) := z(n,m) — z(n,m) (the innovation) is not 0, where z(n,m) is 
the orthogonal projection of z(n,m) on the closed subspace generated by 
{z(k, I) : k < n,l G Z}. In other words, if a deterministic field has non-zero 
column-to-column innovations it cannot be perfectly linearly predicted from 
previous columns. 

When z(n, m) is the deterministic component of the decomposition of 
a regular random field with respect to a NSHP total-ordering, the vertical 
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Figure 2: RNSHP support. 



evanescent component z e (n, m) is the orthogonal projection of z(n, m) on the 
closed subspace generated by the (orthogonal) column-to-column innovations 
{I(k, I) : k < n,l G Z}. Thus, an evanescent field spans a Hilbert space iden- 
tical to the one spanned by column-to-column innovations. In other words, 
the evanescent field is a component of the deterministic field which represents 
column-to-column innovations. Horizontal column-to-column (row-to-row) 
innovations and evanescent components are similarly defined. 

Evanescent processes were first introduced in [7] (on M). In Korezli- 
oglu and Loubaton [TTJ, "horizontal" and "vertical" total-orders and the 
corresponding horizontally and vertically evanescent components of a homo- 
geneous random field on Z 2 are defined. In Kallianpur [9], as well as in 
Chiang [I], similar techniques are employed to obtain four-fold orthogonal 
decompositions of regular (non-deterministic) homogeneous random fields. 
In Francos et. al. [2] this decomposition of random fields on Z 2 was further 
extended. This is done by considering all the Rational Non-Symmetrical 
Half Plane (RNSHP) linear orders, each inducing a different partitioning of 
the two-dimensional lattice into two sets by a broken straight line of ratio- 
nal slope. Intuitively, the usual lexicographic order is not the only possible 
order definition of the 2-D lattice. Each RNSHP linear order is induced by 
a "rotation" of the usual lexicographic order, such that the resulting non- 
symmetrical half-plane is delimited by a broken straight line with rational 
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slope, and which leads to a different linear order definition. Consequently, 
terms such as "past" and "future" are redefined with respect to a specific 
RNSHP linear order (see, for example, Figure [2]). 

More specifically, each Rational Non Symmetrical Half Plane is defined 
in terms of two co-prime integers (a, b), such that the past P a ^ is defined by 

Pa,b = {{n, m) G Z 2 : na + mb < , or na + mb = and m < 0}. (1) 

Then P = P a ^ satisfies 

0') Pn(-P) = {0}, (it) PU(-P) = Z 2 , (Hi) P+P C P (usual addition). 

By (i)-(iii), P induces on Z 2 a linear order, which is defined by (k, I) ^ (n, m) 
if and only if (k — n,l — m) G P. 

Clearly, there are countably many such linear orders. Each such order 
induces a different definition of the term "column", and correspondingly 
different definitions of column-to-column innovations and evanescent field. 

The Wold decomposition of a regular random field into purely-indeterministic 
and deterministic components is invariant to the choice of a RNSHP order. 
The decomposition in [2J further asserts that we can represent the deter- 
ministic component of the field as a mutually orthogonal sum of a "half- 
plane deterministic" field and a countable number of evanescent fields. The 
half-plane deterministic field has no innovations, nor column-to-column in- 
novations, with respect to any RNSHP linear order. On the other hand, 
each of the evanescent fields can be revealed only by using the corresponding 
RNSHP linear order, i.e., with respect to specific definitions of "columns" 
and column-to-column innovations. This decomposition yields a correspond- 
ing spectral decomposition, i.e., we can decompose the spectral measure of 
the deterministic part into a countable sum of mutually singular spectral 
measures, such that the spectral measure of each evanescent component is 
concentrated on a line with a rational slope. 

Based on these results, a parametric model of the homogeneous random 
field was derived in [2J. The purely-indeterministic component of the field 
is modeled by a white innovations driven 2-D moving average process with 
respect to some RNSHP linear order. This component contributes the abso- 
lutely continuous part of the spectral measure of the regular field. One of the 
components of the half-plane deterministic component that is often found in 
practical applications is the 2-D harmonic random field which is the sum of a 
countable number of exponential components, each having a constant spatial 
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frequency and random amplitude. This component contributes the 2-D delta 
functions in the spectral domain. The number of evanescent components of 
the regular field is countable. The model of the evanescent field with respect 
to specific order is presented bellow: 

Let (a, b) be a pair of co-prime integers (a > 0) which defines a specific 
RNSHP linear order according to (jl]). Then, the model of the evanescent 
field which corresponds to this order is 

I(a,b) , v 

e (a,b) ( n i m ) — ( na + mo ) exp ( joo'f^ (nc + md) J , (2) 

i=i ^ ' 

where c and d are co-prime integers satisfying ad — be = 1. For the case 
where (a, b) = (0,1) we have (c,d) = (1,0), and for (a, b) = (1,0) we have 
(c, d) = (0, 1). We further note that in this notation na + mb is the "column" 
index and nc+md defines a "row". The modulating process {s\ a ' b \na + mb)} 
is a 1-D purely-indeterministic, complex valued processes, and oof 1 is a mod- 
ulation frequency. Thus, e( a ,b) has no innovations, with respect to "rows", 
and has non-zero column-to-column innovation (expressed by the modulat- 
ing process s[ a ' b ^) with respect to its "columns". I( a ,b) denotes the number 
of different evanescent components that correspond to the same RNSHP de- 
fined by (a, b). The different components are such that their 1-D modulating 
processes {s^'^} and {s^ a '^}, are mutually orthogonal and their modulation 

frequencies are different ujf 1 ^ ^ ajj' for all 1 < % ^ j < I( a ^ . 

The "spectral density function" of each evanescent field has the form of 
a sum of 1-D delta functions which are supported on lines of rational slope 
in the 2-D spectral domain. The amplitude of each of these delta functions 
is determined by the spectral density of the 1-D modulating process. Since 
the spectral density of the modulating process can rapidly decay to zero, 
so will the "spectral density" of the evanescent field, and hence the name 
"evanescent" . 

1.2 Practical Applications 

Besides its fundamental theoretical importance, the Wold decomposition of 
a regular random field has various applications in image processing and wave 
propagation problems. For example, the parametric model that results from 
these orthogonal decompositions, naturally arises as the physical model in 
problems of texture modeling, estimation and synthesis [3]. 
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Another application is space-time adaptive processing of airborne radar 
data [4J. Space-time adaptive processing (STAP) is an increasingly popular 
radar signal processing technique for detecting slow-moving targets. The 
space dimension arises from the use of array of multiple antenna elements 
and the time dimension arises from the use of coherent train of radar pulses. 
The power of STAP comes from the joint processing along the space and 
time dimensions. Comprehensive analysis of STAP appears in [TOj [T2]. 

In [1] it is shown that the same parametric model that results from the 2-D 
Wold-like orthogonal decomposition naturally arises as the physical model in 
the problem of space-time processing of airborne radar data. This correspon- 
dence is exploited to derive computationally efficient detection algorithms. 
More specifically, the target signal is modeled as a random amplitude complex 
exponential where the exponential is defined by a space-time steering vector 
that has the target's angle and Doppler. Thus, in the space-time domain the 
target contribution is the half-plan deterministic component of the observed 
field. The sum of the white noise field due to the internally generated receiver 
amplifier noise, and the sky noise contribution, is the purely-indeterministic 
component of the space-time field decomposition. 

The presence of a jammer (a foe interference source, transmitting high 
power noise aimed at "blinding" the radar system) results in a barrage of 
noise localized in angle and uniformly distributed over all Doppler frequencies 
(since the transmitted noise is white). Hence, in the space-time domain each 
jammer is modeled as an evanescent component with (a, b) = (0, 1) such that 
its 1-D modulating process s[° (m) is the random process of the jammer 
amplitudes. The jammer samples from different pulses are uncorrelated. In 
the angle-Doppler domain each jammer contributes a 1-D delta function, 
parallel to the Doppler axis and located at a specific angle ujf 1 ' 1 ' using the 
notation of (12]). 

The ground clutter results in an additional evanescent component of the 
observed 2-D space-time field. The aircraft platform motion produces a very 
special structure of the clutter due to the dependence of the Doppler fre- 
quency on angle. The clutter's echo from a single ground patch has a Doppler 
frequency that linearly depends on its aspect with respect to the platform. As 
the platform moves, identical clutter observations are repeated by different 
antenna elements on different pulses, which defines a specific linear locus in 
the angle-Doppler domain, commonly referred as the "clutter ridge" . Thus, 
the clutter ridge, which represents clutter from all angles, is supported on 
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a diagonal line (that generally wraps around) in the angle-Doppler domain. 
Due to the physical properties of the problem the different components of 
the field are assumed to be mutually orthogonal. In the specific applica- 
tion of airborne radar, the evanescent components (the clutter, and jamming 
signals) are considered unknown interferences. 

Although the data collected by STAP radars for different ranges can be 
viewed as a sequence of finite-sample realizations from a homogenous field, 
its is technically more convenient to represent each of the observations in a 
vector form and to statistically analyze them as multivariate vectors. Thus, 
if one uses a STAP system with N antenna elements and M pulses, the 
observed N x M STAP signal is treated as NM x 1 multivariate random 
variable. These vectors are commonly called "snapshots". 

The STAP processor goal is to solve a detection problem, i.e., to estab- 
lish whether a hypothetical target is present or not. It adaptively weights 
the available data in order to achieve high gain at the target's angle and 
Doppler and maximal mitigation along both the jamming and clutter lines. 
The adaptive weight vector is computed from the inverse of the interference- 
plus- noise covariance matrix. [TOl IT2]. It is shown in [5] that the dominant 
eigenvectors of the space-time covariance matrix contain all the information 
required to mitigate the interference. Thus, the weight vector is constrained 
to be in the subspace orthogonal to the dominant eigenvectors. Because 
the interference-plus-noise covariance matrix is unknown a priori, it is typi- 
cally estimated using sample covariances obtained from averaging over a few 
range gates. This is the known as the fully adaptive STAP approach. The 
major drawback of this approach is its high computational complexity. The 
final detection of a target is performed by applying either Constant-False- 
Alarm-Rate (CFAR) detector, or Adaptive-Matched-Filter (AMF) detector, 
or Generalized-Likelihood-Ratio (GLR) detector. Usually the detector is em- 
bedded into the weight computations. 

Fortunately, both the clutter and the jammers have low-rank covariance 
matrices. The clutter covariance matrix has a low rank due to the movement 
of the platform, as discussed above. The jammer covariance matrix has low 
rank since the jamming signal is spatially correlated between all antennas 
at each pulse. The low-rank structure of the interference covariance matrix 
may be exploited to achieve significant reduction in the adaptive problem 
dimensionality with little or no sacrifice in performance relative to the fully 
adaptive case. These methods are referred to as partially adaptive STAP. 

Partially adaptive STAP methods require knowledge of the rank of the 
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interference covariance matrix. However, it is a priori unknown, and unfor- 
tunately cannot be easily estimated from the sample covariance matrix due 
to the existence of a noise component which has a full rank covariance ma- 
trix. Hence, the problem of estimating the rank of the interference covariance 
matrix is critical in the implementation of many STAP algorithms. 

In this work we consider the problem of determining the rank of the 
covariance matrix of a vectorized finite dimension sample from an evanescent 
random field. By using the evanescent field parametric model to model the 
interferences in the STAP problem, we considerably simplify the solution to 
the problem of estimating the rank of the low-rank interference covariance 
matrix. In fact, it turns out that in this parametric framework the well- 
known Brennan rule p2] for the rank of the clutter covariance matrix, as 
well as the rank computation of the jammer, become special cases of the 
general result proved here. Hence, the provided derivation opens the way for 
new, computationally attractive, methods in parametric and non-parametric 
estimation of two-dimensional random fields, with immediate applications in 
partially adaptive space-time adaptive processing of airborne radar data. 

The rest of the paper is organized as follows: In Section[2]we formulate the 
problem we aim to solve. A formula for the rank of the covariance matrix of a 
complex valued evanescent random field is derived in Section [3j In Section [4] 
we extend the obtained result to the case of a real valued evanescent random 
field. Finally, in Section [5] we provide our conclusions. 

The following notation is used throughout. Boldface upper case letters 
denote matrices, boldface lower case letters denote column vectors, and stan- 
dard lower case letters denote scalars. The superscripts (-) T and (-) H denote 
the transpose and Hermitian transpose operators, respectively. By I we de- 
note the identity matrix and by a matrix of zeros. The symbol © denotes 
an element by element product of the vectors (Hadamard product). Given 
a scalar function /(•) and a column vector v, we denote by /(v) a column 
vector consisting of the values of function /(•) evaluated for all the elements 
of v. Finally diag{\) denotes a square diagonal matrix with the elements of 
v on its main diagonal. 
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2 Finite Sample of an Evanescent Random 
Field: Definitions and Problem Formula- 
tion 

Let denote the set of all possible pairs of different co-prime integers (a, b), 
a > 0, where each pair defines a RNSHP order on the 2-D lattice. Although in 
the case of an infinite 2-D lattice the number of different RNSHP definitions is 
infinitely countable, in the finite sample case only a finite number of different 
linear orders can be defined. Moreover, in practical applications the number 
of different evanescent components for each order definition is finite as well. 
Therefore, we assume throughout this paper that \0\ and I( a ,b) are finite 
integers. 

Let Q = Y^r a b)eo h a > b )- Each of the Q evanescent components is uniquely 

defined by the triple (a,b,u>l a ' b ^) where (a,b) G O and uo'f^ is the fre- 
quency parameter. Let us denote the set of all possible triples by Oq = 
{(ai,bi,ux), (ci2, 62, ^2), ■ ■ ■ , ( a Q, bQ, ojq)}. All triples are unique: they either 
have different support parameters (a^&j) 7^ (a,j,bj), or in case (a,,^) = 
(dj, bj) they have different frequencies such that uJi 7^ ujj. 

Finally, adapting ^ to the finite sample case we have that 

Q 

e(n,m) = J^e ff (n,m), (3) 

q=l 

where 

e g (n, m) = s q {na q + mb q ) exp (^juj q {nc q + md q )\ (4) 

such that (a q ,b q ,u q ) G Oq and {s q },c q ,d q are defined as above. 

We note that since the spectral measure of {e q (n, m)} is concentrated on 
a line (that may wrap around) whose slope is determined by a q and b q , we 
interchangeably refer to a q and b q as either the spectral support parameters 
of {e q } or as the RNSHP slope parameters. 

Let {e(n, m) : (n, m) G D} where D = {(n, m) 6Z 2 :fl<n<iV- 1, 0< 
m < M — 1} be the observed finite sample of the random field (|3|. Let e 
denote an NM x 1 vector form representation of this finite sample: 

e = [e(0,0),...,e(0,M-l),e(l,0),...,e(l,M-l),..., 

...,e(iV-l,0),...,e(iV-l,M-l)] r . (5) 
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This is a multivariate representation of a finite sample of an evanescent ran- 
dom field. 

Let r denote the NM x NM covariance matrix of the evanescent vector 



E 



e e 



(6) 



Due to the special structure of the evanescent field, many of the elements 
of e are linearly dependent, and therefore T is low-rank. This property 
is easily demonstrated by considering a single evanescent component that 
corresponds to the vertical order (a, b) = (0, 1) (single jammer source using 
the STAP nomenclature), with some arbitrary modulation frequency u and 
modulating process s(m). In that case 



s{m) exp 



(7) 



The vector form representation of the finite sample of this evanescent field is 



,(o,i) 



[s(0), s(l), . . . , s(M - 1), s(0) exp(j^), . . . , s{M - 1) exp(^), . . . 
s(0) exp(ja;(iV - 1)) . . . , s{M - 1) exp(Ju{N - l))f . (8) 



Since the modulating frequency u is a deterministic constant, it is obvious 
that e^ 0,1 ) is comprised of only M independent random variables. Therefore, 



the rank of r (0ll) = E 



,(o,i) ( e (o,m H 



is also M. 



The aim of this paper is to derive an expression for the rank of the low- 
rank covariance matrix T of the evanescent vector e, in the general case Q). 



3 The Rank of the Covariance Matrix of an 
Evanescent Field 

In this section we derive an expression for the rank of the covariance matrix 
r. In order to do so, we have to find and quantify the linear dependen- 
cies between the samples of e. Unfortunately, for arbitrary spectral support 
parameters and multiple evanescent components, this task involves tedious 
calculations. The results of the entire analysis in this section can be summa- 
rized by the following theorem: 
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Theorem 1. Let e be a vector-form representation of a finite sample from 
a sum of evanescent random fields, given by Then, the rank of its 

covariance matrix, T, is given by 



rank(T) = min NM, 



Nj2\a q \+Mj2\b g \-J2\ a o\Y,\ b ' 

q=l g=l q=l q=l 



<-l\ 



(9) 



Even though the evaluations in the next subsections are technical in na- 
ture, the obtained result is surprisingly interesting. Hence, before addressing 
the proof itself let us make some comments. From Theorem 1, it is clear 
that the rank of the covariance matrix of a finite sample from an evanescent 
random field is completely determined by the spectral support parameters 
(a q , b q ) of the different evanescent components, while it is independent of 
the other parameters of the evanescent fields, such as the parameters of the 
modulating processes, {s q }, or the modulation frequencies, u q . 

Moreover, one can easily observe that the well-known Brennan rule for 
the rank of the low-rank clutter covariance matrix in the STAP framework, 
[T2] as well as the rank of the covariance matrix of the jamming signals are 
special cases of this theorem. The Brennan rule states that the rank of the 
clutter covariance T dut is given by: 

rank(T clut ) = [N + M/3 - (3\ (10) 

where (3 is the slope of the clutter ridge orientation in the angle-Doppler 
domain, and [J denotes rounding to the nearest integer. It is easy to see 
that this formula is a special case of the above theorem when only a single 
evanescent field is observed, and its spectral support parameters are (a, b) = 
The rank of the jamming covariance matrix is, [12]: 

rank(T jam ) = MJ (11) 

where J is a number of sources. Since the spectral support of a single jammer 
in the angle-Doppler domain is a line parallel to the Doppler axis, and since 
all jammers are mutually orthogonal, they can be modeled as J vertical 
evanescent components with spectral support parameters (a, b) = (0, 1), such 
that the rank of the covariance matrix of each individual jammer is M as in 
the above example. 
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3.1 Rank derivations 

In this subsection we prove Theorem 1. The derivation provides an insight 
into the structure of the covariance matrix, and explains the nature of its low- 
rank. Moreover, we explicitly show how columns of the covariance matrix 
that can be represented as linear combinations of other columns are formed, 



which yields its 
Rewriting (3 



ow-rank. 

in a vector form we have e = Yl^=i e 9> where 



e q = [e,(0,0),...,e 9 (0,M-l),e g (l,0),...,e,(l,Af-l), 

..., e<? (iV-l,0),...,e g (iV-l,M-i)] T . (12) 

Let 

= S q( b q)i • • • ' S l(( M - l ) b q)l S gK)> S q( a q + & <?)> ■■■■> S q( a q + ( M ~ l ) b q)' 

. . . , s q {{N - l)a q ), s q ({N - l)a q + b q ), . . . , s q ((N - l)a q + (M — l)b q )f (13) 

be the vector whose elements are the observed samples from the 1-D modu- 
lating process {s q }. Define 



Vq = [ 0, d q , . . . , (M - l)d q , c q , c q + d q ,..., c q + (M - l)d, 



(?;•••; 

(N-l)c q , (N-l)c q + d q , (N-l)c q + (M-l)d q f. (14) 

Let 

D g = diag ^exp(-ju;gV (? )^ . (15) 
be an NM x NM diagonal matrix. Thus, using (|4]), we have that 

e, = Df £ q . (16) 

Let s q be a (N — l)|a 3 | + (M — l)\b q \ + 1 dimensional column vector of 
consecutive samples of the 1-D modulating process {s q }. For the case in 
which a q > and b q < 0, s q is defined as 

s q = [ Sq ((M-l)b q ),...,...,s q ((N-l)a q )f , (17) 

while for the case in which a q > and b q > 0, s q is defined as 

s q = [s q (0),...,... jSq ((N-l)a q + (M-l)b q )f . (18) 
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Thus for any (a q ,b q ) we have that 



and 



q °<3 



V H A T 



where A q is a real-valued [(N — l)\a q 



q ~-q"q 

+ (M 



(19) 



(20) 



l)\b q \ + 1] x NM rectangular 



matrix where each of its columns has a single element whose value is "1", 
while all the others are zero. Thus, each column of A q "chooses" the single 
element from the vector s q that contributes to the corresponding element of 
the vector £ 9 . Due to boundary effects, resulting from the finiteness of the 
observation, not all of the elements of the vector contribute to the vector 
£ , unless \a„\ < 1 or \b q \ < 1. Hence some rows of the matrix A q may contain 



only zeros. On the other hand, whenever na q + mb q 



ka q + 



for some 



integers n,m,k,£ such that < n, k < N — 1 and < m,£ < M — 1, the same 
sample from the modulating process {s q } is duplicated in the elements of £ q . 
Therefore, the number of distinct columns in A q is equal to the number of 
elements of s q that appear in £ , i.e., the number of distinct samples from 
the random process {s q } that are found in an observed finite sample of an 
evanescent field of dimensions N x M. The matrix A q depends only on (a q , b q ) 
and is independent of the modulation frequency u q or the modulating process 

The rank of covariance matrix T is strongly related to the number of 
distinct samples from the random processes {s q } for all 1 < q < Q which 
can be found in the evanescent vector e. Therefore, the rank of T is tightly 
related to the ranks of the matrices A q , 1 < q < Q. 

Let R g be the covariance matrix of the vector s q i.e., 



H q = E 



(s,) 



(21) 



The matrix H q is full rank positive definite since the process {s q } is purely- 
indeterministic. Since the evanescent components {e q } are mutually orthog- 
onal we conclude that T, the covariance matrix of e, has the form 



T = E 



(e) 



H 



9=1 



(22) 
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where T q is the covariance matrix of e q . Using (20) and (21) we find that 



T q = E 



e q \ e qJ 



Df A^R g A g D g . (23) 



Finally, 

Q 

r = ^DfAjR g A g D g . (24) 

9=1 

One can rewrite the above expression in a block-matrix form 

T = CJRCq (25) 

where 

C =[DfA?...D£A5]*, (26) 

and 

R = dia^([Ri...RQ]). (27) 

is a block-diagonal matrix with the matrices R g , 1 < q < Q on its diagonal, 
and zeros elsewhere. 

Since the covariance matrices R g are full rank positive definite, the block- 
diagonal matrix R is full rank positive definite as well. Hence by observation 
7.1.6 [S] we have 

rank{T) = rank (Cq) . (28) 

The matrix Cq has exactly NM columns, such that each one of its 
columns corresponds to an entry in the evanescent vector e, or similarly, 
each one of its columns corresponds to a point on the original N x M lattice 
D = {(n,m) E Z 2 : < n < N - I, < m < M - I}. More specifically, 
the n(M — 1) + m column of Cq corresponds to the n(M — 1) + m element 
of e which represents the evanescent field sample at the (n, m) lattice point. 
(Note that we enumerate the columns starting from zero). In the following 
we will adopt the abbreviation [n, m] for indexing the n(M — 1) + m column 
of a matrix. 

To gain more understanding on the structure of Cq let us examine the 
different matrices Cq is comprised of. We begin with A q for some 1 < q < Q: 
By construction (see (20) and the following explanation) all columns of A q 



are unit vectors, where the single "1" entry in each column chooses the single 
element from the vector s q that contributes to e(n, m) - the evanescent field 
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sample at (n, m) . The single non-zero entry in the [n, m] column of A q is 
located in the fc'th row where na q + mb q = k (we allow negative indexed rows 
in the case where b q < 0). For example if a q > and b q > 0, the matrix A q 
is given by 



[0,0] 
1 



[N 



k 



(N - IK + {M — 1)6, 



1,M- 
■■ 



1] 











(29) 



Let (n* 



m 



be a solution to the linear Diophantine equation na q + mb, 



g 



ta 



k. Then, the equation is also satisfied by n = n* + tb q and m 
where t is an arbitrary integer. Since (a q , b q ) are coprime integers these are 
the only possible solutions. It means that as soon as (n + tb q , m — ta q ) G D, 
the corresponding [n + tb q , m — ta q ] column of A q will be equal to its [n, m] 
column. To find the rank of A, we have to evaluate the number of linearly 
independent columns, i.e., the number of distinct elements of s q (na q + mb q ) 
(n, m) G D which contribute to e. 

Since D g is a diagonal matrix, the structure of A g D g is similar to the 
structure of A q with the only difference being that instead "1" in each col- 
umn, we have the appropriate exponential coefficient. Therefore, each col- 
umn of the matrix Cq has exactly Q non-zero elements. 

Next, let us concatenate the matrices A p and A q , where 1 < p ^ q < Q 
and examine the structure of resulting matrix 



'pq 



[A p A q ] . 



(30) 



As before, let us consider the structure of some arbitrary [n, m] column of 
this matrix. It has two non-zero entries: On the k p row of A p and on the k q 
row of A q , where (n,m) satisfies 



na p + mbp 



k 



(31) 



and 



na„ + mbg 



kq. 



(32) 



Next, we note that the pair (n + tb p , m — ta p ) satisfies the linear Diophan- 



tine equation (31) for any integer t. Therefore, for t p such that (n + t p b p ,m- 
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tpdp) G D, the [n + t p b p , m — t p a p } column of A p has a "1" entry, at the same 
row as the [n, m] column. However, (n + t p b p , m — t p a p ) also satisfies the 
linear Diophantine equation 



(n + t p b p )a q + (m — t p a p )b q = I, 



(33) 



Hence, the [n + t p b p , m — t p a p ] column of A q has a "1" entry on its i q row. 

Similarly, since (n + tb q ,m — ta q ) satisfies the linear Diophantine equation 
(32) for any integer t, for t q such that (n + t q b q , m — t q a q ) G D we have that 
the [n + t q bq, m — t q a q ] column of A q has "1" at the same row as the [n, m] 
column. Since, 

(n + t q b q )a p + (m - t q a q )b p = £ p , (34) 



the [n + t q b q ,m — t q a q ] column of A p has "1" on its £ p row. Moreover, one 
can observe that for a pair of integers (t p , t q ) such that {n + t p b p + t q b q , m — 
t p a p — t q a q ) G D, the pair {n + t p b p + t q b q , m — t p a p — t q a q ) simultaneously 



satisfies (33) and (34): 



(n + t p b p + t q b q )a p + (m - t p a p - t q a q )b p 
(n + tpbp + t q b q )a q + (m - t p a p - t q a q )b q 



(35) 



Therefore the [n + t p b p + t q b q , m — t p a p — t q a q ] column of of A p has "1" on 
its t v row, and the same column of A q has "1" on its i q row. 

Finally, we can represent the [n, m] column of C pq by a linear combination 
of its other columns: 



n, m\ 



[n + tpbp, m - t p a p ] + [n + t q b q , m - t q a q ] - 
[n + tpbp + tqbq, m - t p a p - t q a q ], 



(36) 



or in a more detailed form by 

1 
1 



+ 



(37) 



Let Tp q ' m ^ be the set of all the integer pairs (t p , t q ) such that (n + tpb p ,m — 
tpdp), {n + t q b q , m — t q a q ), {n + t p b p + t q b q , m — t p a p — t q a q ) G D. Clearly, the 
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set Tpq' mS) is non-empty since (0,0) G T pq ' m \ and it corresponds to a trivial 
representation of the [n,m] column by itself. If \T pq ' m ^\ > 1 then the [n, m\ 
column has non-trivial linear representation by other columns. 

Recall however, that the matrix Cq is comprised of blocks where each 
block is of the form A g D g . Consider next the concatenation of two such 
blocks ApDp and A q D q , 

C pg = [Df AJ Df A T q } H . (38) 

Keeping in mind that by definition a p d p — b p c p = 1 and a q d q — b q c q — 1, it is 
easy to check that the replacement of the "1" in the columns of C pq by expo- 
nentials as in (38) will only affect the coefficients of the linear combination. 
Indeed, the linear combination of columns in this case has the form 

[n, m] = [n + t p b p , m - t p a p ] exp(ju p t p ) + [n + t q b q , m - t q a q ] exp(ju q t q ) 
- [n + t p b p + t q b q ,m-t p a p -t q a q }exp(j[uj p t p + uj q t q ]). (39) 
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One may also notice that if (a p , b p ) = (a q , b q ) we have that k p = k. 
However, since in this case u p ^ u q the linear combination in (39) is still 
valid and non-trivial. 

It is clear that the linear dependencies of columns of Cq are governed 
by the same simple laws: Let Tq" 1 ' be a set of all Q-tuples of integers 
(ti, . . . , to) defined as follows: For any 1 < q < Q, let (ii, be a set of 

q indices, such that 1 < < Q for all 1 < k < q, and 

(n + t h b h + . . . + t iq b iq ,m- t h a h - ... - t iq a iq ) E D. 

Clearly, the set Tq 1 '" 1 ^ is non-empty since (0, . . . , 0) G Tq ,171 \ Let (n, m) G D 
be an arbitrary lattice point and let [n, m] be its corresponding column in 
Cq. Then, the [n,m] column can be represented by the linear combination 



Q Q-q+l Q 

^(-1) 9_1 ^ ... ^ [ n + t i 1 K + ••• + ti g bi g ,m-t il a h - ...-t ig a,i 

q=l il = l jq=jq_l + l 



v 

q sums 



exp {]{ui 1 t h + .. . + uj iq t iq }) , (40) 



where (ti, . . . ,tq) G Tq™'" 1 ^. The details of this derivation are presented in 
Appendix A. 
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Following the foregoing analysis of the linear dependencies between the 
columns of Cq, we next count its linearly independent columns in order to 
derive the rank of Cq. Let us first count the number of independent columns 
of AqDq. As mentioned earlier, this number is equal the number of distinct 
samples from s q (na q + mb q ), (n,m) G D that contribute to e. In other 
words, this is the number of different indices k, such that na q + mb q = k 
where (n, m) G D, and it can be easily calculated based on the dimensions 
of D (see Figure [3] for an illustrative example). Indeed, a new sample from 
the random process {s q } may be introduced only on the first a q rows (since 
a q > 0) and the last (first) \b q \ columns (last if b q > and first if b q < 0) 
of the observed finite dimensional field, while the rest of the field is filled by 
replicas of these samples. We thus count Na q distinct samples in the first a q 
rows and M|6 9 | distinct samples in the first (last) \b q \ columns. However, on 
the intersection of these rows and columns \a q b q \ samples are counted twice. 
Finally, the total number of distinct samples from the random process {s q } 
that are found in an observed field of dimensions N x M (which is equal to 
the rank of A q and the rank of A q T) q ) is given by 

r q = Na q + M\b q \ - \a q b q \. (41) 

Similarly, it can be shown that the number of linearly independent columns 
of ApDp is Na p + M \ b p \ — \ a p b p \ . Let us next count the number of linearly 
independent columns of C pq . 

Since r p columns of A p D p are linearly independent, the same columns of 
the concatenated matrix C pq are linearly independent as well. The remaining 
NM — r p = (N — \b p \)(M — \a p \) columns may be considered to correspond 
to an (N — \b p \) x (M— \a p \) rectangular sub-lattice Di = {(n, m) G Z 2 : < 
n < N — 1 — | ftp |, | a p | < m < M — 1} (see Figure [4] as an example), which 
is a subset of the original rectangular lattice (or similarly, one can define 
D 1 = {{n,m) G Z 2 : \b p \ < n < N - 1, |a p | < m < M - 1} which doesn't 
change the reasoning of our arguments and only depends on a sign of b p ). 
Repeating the same arguments as those made above, one can show that the 
number of distinct samples from the random process {s q } that are found in 
a sub-lattice Di is 

r~ = (N — \b p \)\a q \ + (M — \a p \)\b q \ - \a q b q \ . (42) 

This is the number of linearly independent columns which can be found in 
C pq in addition to the first r p columns. 
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Figure 3: N = M = 15, a = 3, b = ±2: (a) The indices k = na + mb of the 
observation on {s( 3 ' 2 )(/c)}; (b) The indices k = na + mb of the observation 
on {s^'^(k)}. The sets of distinct indices of {s^ 3 ^(k)} and {s( 3 ~ 2 )(/c)} are 
marked in yellow. Every other sample in the field is identical to some sample 
in the yellow area. 
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Figure 4: N = M = 15, a x = 3, 6 a = 2, a 2 = 2,b 2 = ±1. (a): The set of 
distinct samples of s^ 3 ' 2 ^ (n, m) (in yellow) and s^ 2 '^(n, m) (in blue), (b): The 
set of distinct samples of s( 3,2 )(n, m) (in yellow) and s^ 2,_1 ^(n, m) (in blue). In 
both cases, every other sample in the field is a linear combination of samples 
in the colored areas. 
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Let D 2 be the set of NM - r p - f q = (N - \b p \ - |6 P |)(M - \a q \ - \a q \) 
lattice points that remain after the removal from D of the r p + r q points 
corresponding to the linearly independent columns of C pq (for simplicity and 
without limiting of the generality of the results, we will discuss the case 
where b p > and b q > 0, as illustrated in Figure 2a (uncolored area)). Thus, 
D 2 = {{n,m) EZ 2 :0<n<N-l - \b p \ - \b q \, \a p \ + \a q \ < m < M - 1}. 
It thus remains to be shown that all columns representing points in D 2 can 
be represented by a linear combination of columns that correspond to points 
in D \ D 2 . 

Since the "width" of D \ D 2 is \b p \ + \b q \ along the n-axis and \a p \ + \a q \ 
along the m-axis (colored areas in Figure 4), for every (n, m) G D 2 we have 
{n + b p ,m — a p ), {n + b q ,m — a q ), {n + b p + b qi m — a p — a q ) G D. Thus, 
(t p ,t q ) = (1,1) G T pq ' m \ and as we have shown above, we can represent 
[n, m] by the linear combination 

[n, m] = [n + b p , m - a p ] exp(juj p ) + [n + b q ,m- a q ] exp(juj q ) 

- [n + b p + b q ,m- a p - a q ) exp(j[u p + u q )). (43) 

Continuing this construction recursively, it is obvious that for each point 
(n, m) G D 2 we can find a pair (t p , t q ) G T pq ' m \ and (t p} t q ) ^ (0, 0) such that 
(n+t p b p ,m-t p a p ), (n+t q b q ,m-t q a q ), (n+t p b p +t q b q , m-t p a p -t q a q ) G D\D 2 . 
On the other hand, for every point (n, m) G D\D 2 one can show that T pq ' m ^ = 
{(0,0)}, i.e., only the trivial linear combination exists. In other words, all 
the random variables indexed on D \ D 2 correspond to linearly independent 
columns. Therefore, the number of linearly independent columns in C pq is 

rank(C pq ) — \D \ D 2 \ = r p + f q 

= N(\a p \ + \a q \) + M(\b p \ + \b q \) - (\a p \ + \a q \)(\b p \ + |6,|). (44) 

The construction described above can be easily extended to the general 
case where we concatenate all the matrices which Cq is comprised of. See 
Figure [5] for an example of a three component case. If one chooses a subset 
of the original lattice, 

Q Q 

D Q = {(n,m) eZ 2 :0<n<N-l Y^K| < m < M - 1}, 

which remains after the removal of 

Q Q Q Q 

Nj2\a q \ H E \ b «\ ( 45 ) 

q=l q=l q=l q=l 
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Figure 5: Sets of distinct samples in the case of three evanescent components. 
N = M = 15, ai = 3, b x = 2, a 2 = 2,b 2 = 1, a 3 = 1, b 3 = 3 



lattice points (similarly to D 2 which remains after the removal of the r p + r q 
lattice points corresponding to the independent columns of C pq ), one can 
repeat the same considerations as above and show that columns of Cq that 
correspond to the lattice points in D \ Dq are the only linearly independent 
columns of Cq. Thus, 

Q Q Q Q 

rank(C Q ) = \D \ D Q \ = Nj^l^l + M J2 IM ~ E l a «l ^ 

q=l g=l q=l q=l 

Finally, Since the rank of T cannot exceed NM (the dimension of the 
covariance matrix), NM is an upper bound on the rank of F. Combining 



this and (46) the rank of T is given by fl9h, which completes the proof. 



4 The Case of a Real Valued Evanescent Field 

In the case where a real valued evanesced field is considered, we have 

e q (n, m) = s q {na q + mb q ) cos{ui q {nc q + md q )) 

+ t q (na q + mbq) sm(u q (nc q + md q )), (47) 
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where the 1-D purely-indeterministic processes {s q }, {s p }, {t q }, {t p } are 
mutually orthogonal for all 1 < p 7^ q < Q, and for all q the processes 
{s q } and {t q } have an identical autocorrelation function. Let t q be defined 



similarly to s„ in (18). Using similar notations as in (20) we have 



e g = ^(Df)A^s g +J(Df)A^t 9 , 



(48) 



where 71 and X denote real and imaginary parts respectively. Finally, since 
the processes {s q } and {t q } are mutually orthogonal and have an identical 
autocorrelation function we find that 



E 



e <? \ e qJ 



H(D? )AjB,A,W(D fl ) + Z(Df )AjH^I(D,) (49) 



where 



R„ 



E 



S q l S d 



E 



t q (tq) 



(50) 



is positive definite since {s q } and {t q } are purely-indeterministic. 

Similarly to the case of a complex valued evanescent field, the covariance 
matrix T is given by 



9=1 



(51) 



The derivation of the rank of the covariance matrix (51) follows exactly 



the same lines as in the previous section, and the next corollary is immediate: 

Corollary 1. Let e be a vector-form representation of a finite sample from 
a sum of real-valued evanescent random fields, given by Q), (J^l) and 
Then, the rank of its covariance matrix, F, is given by 



rank(T) = min NM, 



q=l q=l q=l q=l 



(52) 



5 Conclusion 

We have considered the problem of evaluating the rank of the covariance 
matrix of a finite sample from an evanescent random field. We have an- 
alytically derived the rank formula and have shown that the rank of the 
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covariance matrix of this finite sample from the evanescent random field is 
completely determined by the evanescent field spectral support parameters 
and is independent of all other parameters of the field. Thus, for example, 
the problem of evaluating the rank of the low-rank covariance matrix of the 
interference in space time adaptive processing (STAP) of radar data may be 
solved as a by-product of estimating only the spectral support parameters 
of the interference components, when we employ a parametric model of the 
STAP data which is based on Wold decomposition, [3]. Thus , this formula 
generalizes a well known result known as the Brennan rule for the rank of 
the clutter covariance matrix in space-time adaptive processing of airborne 
radar data. The derived rank formula may be employed in a wide range of 
applications in radar signal processing as well as in other areas of signal and 
image processing. 

6 Appendix A 

To derive ( |40| ) let us choose an arbitrary lattice point (n,m) E D, and hence 
a corresponding column [n, m\. Exactly as in the two component case, this 
column is associated with the random variable e(n, m). In fact we are looking 
for a representation of this random variable by a linear combination of other 
random variables indexed on D. 

The first term in the desired linear combination is a sum of Q columns 

Q 

^J[n + Ubi, m - Uai] exp(juiU). (53) 
i=i 

This sum creates a new column vector. Similarly to the two component 
case, this vector is composed of the Q non-zero elements of the [n, m] column 
(similarly to the elements in rows k p and k q of C pq ), but in addition it includes 
the undesired elements (similar to the elements in rows i v and £ q ). The 
total number of contributed undesired elements is Q(Q — 1). Each two pairs 
(aj,6j) and (dj,bj), i ^ j, contribute a pair of undesired elements (one in 
AjDj and one in Aj~Dj), which can be eliminated by subtraction of the 
[n+tibi+tjbj,m—tiai—tjaj] column multiplied by the appropriate exponential 
coefficient, since 

(n + tibi)cij + (m — tidi)bj — (n + tjbi + tjbj)dj + (m — tjOj — tjdj)bj 

{n + tjbj)di + (m — tjdj)bi = (n + tibi + tjbj)di + (m — Udi — tjdj)bi (54) 
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(See also (35) for the equivalent scenario in the two component case). 

To eliminate all these undesired elements we subtract from the vector in 
(53) the sum of (^) such columns (half the number of contributed undesired 



elements), and the result is 

Q 

^[n + Ubi, m - Uai] exp(juiti) - (55) 
i=i 

Q-l Q 

/ j 2J i n + Uh + tjbj, m - tidi - tjdj] exp (j[u;iU + uJjtj)) 

i=l j=i+X 

However, the resulting column now contains a new kind of undesired ele- 
ments. Substraction of the [n+tibi+tjbj,m—Uai—tjaj] column has eliminated 
2 undesired elements in AjDj and AjDj, but at the same time has created 
a new undesired element in every A^D^, such that i 7^ j 7^ k. A total of 
Q — 2 new undesired elements have been created. Clearly, these elements are 
negative and their total number is (®) (Q — 2) =3 (^) . To eliminate the con- 
tribution of these elements we add the [n+tibi+tjbj+t k b k , m—tiai—tjaj—ti-ak) 
column with an appropriate exponential coefficient which eliminates the un- 
desired element from A^Dfe. At the same time this action is also canceling the 
undesired result of subtracting [n + tibi + tkbk.m — tiai — tkak] that appears in 
AjDj, and the undesired result of subtracting [n + tjbj + tkbk,m — tjaj —tkdk] 
in AjDj. This is because it can be easy verified that indeed 

in + tibi + tjbj)dk + (m — t^di — tjdj)bk = 

(n + + tjbj + t k b k )a k + (m - - tjdj - t k a k )b k , 
(n + Ubi + t k b k )aj + (m - - t k a k )bj = 

(n + t^i + tjbj + tkbkjdj + (m — tjOj — tjdj — tkdk)bj, 
(n + tjbj + t k b k )ai + (m - tjdj - t k a k )bi = 

(n + Ubi + tjbj + t k b k )ai + (m - Uat - tjaj - t k a k )bi. (56) 



To eliminate all undesired elements we add to the vector in (55), (3) such 
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columns, and the result is 

Q 

^2\n + Uh, m - Uai] exp(jujiti) - 
i=i 

Q-l Q 

^ [ra + Ubi + tjbj, m - Uai - tjdj] exp (j[uJiU + ujjtj]) + 

i=l j=i+l 
Q-2 Q-l Q 

^2 Y Y ^ + tibi + t ^ > 3 + ^ &fe ' m ~~ ti<Xi ~ t i <1 3 ~ tk<lk ^ ' 
i=l j=i+l k=j+l 

■ exp (j [uiU + ujjtj + uj k t k }) (57) 

The last action canceled (^)(Q — 2) = 3(^) undesired elements and cre- 
ated once again new (^) (Q — 3) =4 (^) undesired elements. We repeat 
this procedure Q times and in each step k, 1 < k < Q, we subtract/add 
(2) columns for canceling k(®) undesired elements created in the previous 
step. Due to this substraction/addition new (^)(<5 — k) = (k + l)( fc ^ 1 ) un- 
desired elements are created. Clearly, when we subtract/add columns 
exactly Q = Q (^) undesired elements are created. These may be canceled by 
substraction of a single vector. By substraction/addition of the last vector, 
[n + t\bi + . . . + tg6Q,m — t\ai — ... — tQdQ] the process terminates, since 
we are canceling the last Q undesired elements and remain with Q elements 
- exactly those of the [n,m] column, i.e., 

Q 

[n, m] = ^2[n + Ubi, m - Uai] exp(juiU) - 
i=i 

Q-l Q 

[n + Uhi + tjbj, m — U^i — tj a j] ex P Ui^iU + ^jtj]) + • • • 

i=l j=i+l 

+ t 1 b 1 + . . . + t Q b Q ,m- Mi - . . . -t Q a Q ] exp (j^U + . . . + u Q t Q ]) 
Q Q-q+l Q 

X^ -1 )^ Yl ■■■ Yl [n + t il bi 1 + ... + ti q bi q ,m-ti 1 a il -...-ti q ai q \- 

q=l h=l i q =i q - 1 +l 

V v ' 

q sums 

•exp(j[cJiA + . . . + u iq t iq \) (5 
Clearly, this linear combination will be meaningful only if (t\, . . . , £q) G 
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T^ ,m> , i.e., (n + t h b h + . . . + t iq b iq ,m - t h a h - ... - t iq a iq ) e D, for any 
1 < Q < Q, and where . . . , i q ) is such that 1 < ik < Q for all 1 < k < q. 
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